Stochastic description of the dynamics of a random-exchange Heisenberg chain 
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We study the diffusion process in a Heisenberg chain with correlated spatial disorder, with a power 
spectrum in the momentum space behaving as k~'^ , using a stochastic description. It establishes a 
direct connection between the fluctuation in the spin-wave density of states and the noise density 
of states. For continuous ranges of the exponent /3, we find super-diffusive and ballistic spin- 
wave motions. Both diffusion exponents predicted by the stochastic procedure agree with the ones 
calculated using the Hamiltonian dynamics. 



I. INTRODUCTION 



In the last decades, a considerable number of dynamical systems has been studied and a great deal of attention has 
been paid to the analysis of its transport properties. In particular, the study of diffusion and transport properties of 
physical systems with short or long-range correlations in the disorder distribution has attracted a renewed interest |l|, 
QiHi^lEIEI^IEIE^iOiQ' For instance, the unexpected high conductance of several doped quasi-one-dimensional 
polymers was explained by assuming pairwise correlations in the disorder distribution P| . Similarly, the suppression of 
Anderson localization was recently confirmed experimentally in semiconductor superlattices with intentional correlated 
disorder Further, it was demonstrated that long-range correlations in site also act towards the delocalization of 
ID quasiparticle states The ID Anderson model with long-range correlated diagonal disorder displays a finite 

phase of extended states in the middle of the band of allowed energies, with two mobility edges separating localized 
and extended states |lClj |. This result was experimentally validated by microwave transmission spectra through a 
single-mode waveguide with inserted correlated scatterers .14] . The above results have motivated the study of other 
model systems that can be mapped onto the Anderson model, such as magnetic ^lE] and harmonic chains 16]. 

In the context of stochastic process, Morgado et al. 0, studied diffusion in systems with long-range time correlation. 
First, they establish a direct connection between the noise density of states, yOn(w), and the diffusive process. Second, 
they conjecture that the dynamics of a Hamiltonian system with space correlated disorder could be described by 
the same formalism if one supposes that the fluctuation in the density of states of the quasi-particle or elementary 
excitation, pf{oj), plays the same role as p„ in the stochastic description. In this work, we present a numerical 
analysis of the validity of this conjecture. We study it for the one-dimensional quantum Heisenberg model exhibiting 
long-range correlations in the random exchange couplings. For continuous ranges of the degree of correlation, this 
system presents super-diffusive and ballistic motions 15]. Here, we provide numerical evidence that the Hamiltonian 
description and the stochastic one can be unified through the referred conjecture, thus confirming early expectations 

Q. 

II. STOCHASTIC AND HAMILTONIAN DESCRIPTIONS 

Let us suppose that the equation of motion for an operator A can be cast in the form 0, 0, 

dA{t) 



dt 



r{t-t')A{t')dt' + h{t), (1) 
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where h(t) is a stochastic noise subject to the conditions {h{t)) = 0, {h{t)A{0)) = 0, and to the fluctuation-dissipation 
theorem |l7l |: 

Ch{t) =< h{t)h{Q) >=< > r{t). (2) 

Here, (...) indicates an ensemble average in thermal equilibrium. In principle, the presence of the kernel T(t) allows 
us to study a large number of correlated processes. For example, in analogy with the usual Langevin's equation, we 
can study the asymptotic behavior of the second moment of the variable, 

a(t) = [ A{s)ds, (3) 







namely 



= (4) 



where if is a constant. In Eq. Q), we have a = 1 for normal diffusion; for sub- and super-diffusion, a < 1 and a > 1, 
respectively. 

The generalized field, h{t), in Eq. Q can be modeled by a thermal bath composed of harmonic oscillators; conse- 
quently, according to Eq. the memory function can be written as 

T{t) — J pn{oj) coa{ujt)duj, (5) 

where p„(a;) is the noise density of states. It has been proved 0] that If the Laplace transform of the memory function 
of a unidimensional system behaves as 

f (z ^ 0) cx z", (6) 

then the diffusion exponent is 

a = v + l. (7) 

In disordered Hamiltonian systems the diffusion process can be studied through direct integration of the equations of 
motion 0,^^. Now then, how can we assure that the two approaches are compatible and lead to the same results? 
For those systems the density of states of the quasi-particle or of the elementary excitation, D(i?), plays the most 
significant role. However, it displays fluctuations, which are intrinsically connected to the diffusive behavior. Here, 
we conjecture that, for the relaxation processes, the fluctuation in the density of states, pp{E), plays the same role 
as the noise density of states in the stochastic process, thus 

p,,{E) ^ pf{E). (8) 

Consequently, if we calculate pf{E) we can obtain the memory function through Eq. (O and the diffusion exponent 
using Eqs. lO and ITJ. Below, we explore these ideas in the framework of the referred magnetic Hamiltonian system 

. . . . . 

The one-dimensional Heisenberg chain with long-range correlated random exchange couplings can be described by 
the Hamiltonian 

N 



H = -J2Ji^i-Si+i, (9) 



1=1 

where S = 1/2, and the exchange integral is deflned by 

N/2 

J; = ^A(fc)l/2cos(7rfcViV-f 0fe). (10) 
k=l 

Here, (j)k are random phases, and A(k) is a power law spectrum in the k space given by 

A(A:)cxfc-''. (11) 



FIG. 1: Fluctuations in the density of states pf{E) as function of E for /? = and P — 1.5. The average is over 900 samples. 
For /3 = 0, the curve suggests a power-law for low energies, while for P — 1.5, a plateau is observed. 



By using a numerical renormalization group technique, it was predicted that this Hamiltonian supports a phase of 
low-energy extended spin waves in the strongly correlated regime (3 > 1. Associated with these non-scattered modes, 
the spread of an initially localized wave packet displays a ballistic behavior in the long time limit i.e. Eq. Q 
with a = 2. In the weakly correlated regime, < /? < 1, a super-diffusive behavior is obtained with a = 1.5. It is 
worth mentioning that the wave packet dynamics was investigated by means of a direct integration of the equations 
of motion using a fourth order Runge-Kutta method 0,^3- ^^^^ <^^se the site spin operator S; is here identified 
with operator A in the stochastic formalism. Therefore if conjecture (jS)) is true the exponent of Eq. {Tj) should match 
the one calculated using the Hamiltonian dynamics |15|. 



III. NUMERICAL ANALYSIS 

Now we describe the numerical method used to obtain the dynamical behavior, which consists in calculating the 
Laplace transform r{z) of the memory fmiction T{t) defined by Eq. However, we first compute the function 
Pf{E), defined by 

Pf{E) T){Ef >c-< r>{E) >l, (12) 

where is an average over chains with different random sets of exchange integrals. To that end, we use Dean's 

numerical method |2lj | to obtain the spin- wave density of states, D(i?), and the corresponding pp{E), shown in Fig. 
1, for /3 — and (3 — 1.5. It is clear from the data that, for low energies, the fluctuations pf{E) in the correlated 
case, (3 = 1.5 in this figure, are smaller than the uncorrelated case (/? = 0). For /3 > 1, it was shown in ref 15^ that 
the low energy states are extended, and that the density of states has a similar behavior as that of the pure chain. 
Therefore, we concluded that the finite value for pf{E) obtained in the correlated case {(3 = 1.5) was not stable in 
the thermodynamic limit. In order to obtain the most probable value of the fluctuations in the correlated case, we 
perform a scaling analysis of the function pf{E) for chains with distinct lengths. 
Hence, we consider the following scaling function 

Q{E,Ni,N2) ^ eMAPF{E,Ni)Ni ~ Pf{E,N2)N2\), (13) 

with iVi > N2- This method has been successfully used to obtain the behavior of a harmonic chain in the thermo- 
dynamic limit |1^. In Fig. |21 we show the scaling as a function of the energy E. If Q{E, Ni, N2) vanishes when we 
increase the chain size, the fluctuations are flnite in the thermodynamic limit. However, if Q{E, Ni, N2) ~ 1, the 
fluctuations pf{E) decrease to zero when we increase the chain size. In order to clarify this figure, we interpolated 
the data using a Bezier interpolation function. In the insets we show the original data without interpolations. We 
use A'2 = 1.0 ■ 10"*, and three distinct values for Ni. Fig. 2(a) displays the results for /3 = 0. From this figure we can 
see that Q{E, Ni, N2) ^ 1 in all energy ranges, except at i? = 0. This indicates the existence of an extended state 
at this energy, which is determinant in the super-diffusive behavior. Figure|2fb) shows Q(E, Ni, N2) for /? = 1.5. In 
the range of low extended states, < E < Ec, Q{E, Ni, N2) — > 1, meaning that pf{E) in the thermodynamic 
limit. Above Ec, Q{E, Ni, N2) ^ 1. The presence of a finite range of energy with extended states, pf{E) = for 
< E < Ec {Ec ~ 4) , is responsible for the ballistic behavior . 
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FIG. 2: Scaling function Q{E, N\, N2) as function of energy E. Here 7V2 = 1.0 • 10'' in all figures, (a) From top to bottom 
iVi = 1.5 ■ 10", 2.0 • 10*, 3.0 • 10*. We find that 6^1, i.e. pf{E) = only for E = 0, otherwise 9 ^ in the thermodynamical 
limit. The inset shows the original data without interpolations for A'^i = 3.0 • 10*. Average performed over 120 chains, (b) 
From top to bottom A'^i — 2.0 ■ 10*, 3.0 • 10*, 4.0 • 10*. In the thermodynamical limit O ^ 0, except in the region < E < Ec 
{Ec ~ 4) where pf{E) and 0^1. Inset: A'^i — 4.0 • 10*. Average performed over 100 chains. 



Now we proceed by numerically integrating pf{E) in order to obtain the limiting behavior of r(z) as 2; 0, 
using Ec = 4. In Fig. |2| we show T{z) as a function of z for (3 — 1.5. Observe the verification of the conjecture: 
V Ki 1.0, implying a » 2.0. We find super-diffusive motion for [3 < 1, and ballistic motion for /? > 1. Superdiffusive 
propagation of a wave packet was also obtained in the one-dimensional Fibonacci and Thue-Morse lattices, which 
present deterministic aperiodic order. Since the numerical method has greater errors for smaller values of we find 
it difficult to obtain the density of states for small E (large times). This sort of error is expected to have an influence 
on small values of z in the Laplace transform. Another numerical problem is due to the finite size of the chain, which 
has important implications for /3 > 1, as seen by the presence of the plateau in Fig.^ 

Finally, Fig. ^ shows a as a function of 13. We select < /? < 2, covering both regimes. The numerical simulations 
are quite compatible with the stepwise behavior. The maximum deviation observed is at /3 = 0.8, with an error of 
about 8%, which results in a = 1.39. 



IV. DISCUSSIONS AND CONCLUSIONS 



The main results reported in this work are the following: first, for the Heisenberg system we have calculated the 
superdiffusive and ballistic exponents with reasonable precision. Second, we have been able to associate a memory 
function to a Hamiltonian, through which we unite two powerful formalisms: a stochastic description and the quantum 
equations of motion. This result may have implications far beyond the specific example outlined in this work. 
Evidently, we need more simulations and experiments connecting the exponents a and v in several diffusive processes 
in order to have a more complete picture of the validity of conjecture ©. It is important to notice that there are 
many phenomena associated with those described here. For example, in the ballistic regime (a = 2), there remains 
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FIG. 3: The Laplace transform of the memory function r(2:) as function of z. The upper curve is a guide for the eye. The lower 
curve is our result for /3 — 1.5, averaged over 900 chains. In the limit z we approach the expected linear behavior. 
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FIG. 4: Numerical data for the diffusion exponent a as function of (3 averaged over 900 chains. The step is the expected 
result O. 

some open issues concerning the use of the fluctuation-dissipation theorem ^251. Since the early experiment of 
Kauzmann [24|. situations have been found where the systems do not thermahze, usually associated with their slow 
dynamics |22l \24. l25l . This slow dynamics is due to the softening of the lower fluctuating modes, here 

identified by the presence of a finite range of energy with extended states: pf{E) = for < E < E^- This effect 
has been used recently to obtain ballistic diffusion in ratchet devices. Another related phenomenon is chaos 
synchronization |30j |. which presents a reduction in the allowed phase space. A large time to achieve equilibrium, or 
the absence of some re gion s of the phase space, leads to the same effect. Some of its important consequences, such as 
violation of ergodicity |22l Isif and of the fluctuation-dissipation theorem , have been observed. 
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